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Are  there  principles  to  guide  eddy  parameterizations? 

Greg  Holloway 

Institute  of  Ocean  Sciences,  Sidney,  B.C.,  Canada 


Abstract.  A “list  as  long  as  your  arm”  has  described  the  range  of  principles  that 
may  be  attempted  for  possible  eddy  parameterization  schemes.  We  feel  uncom- 
fortably an  attitude  of  “try  and  see”  whether  any  particular  principle  “works”  in 
any  particular  application.  This  note  summarises  a discussion  which  followed 
the  body  of  presentations  at  ‘Aha.  We  ask  if  principles  from  general  physics, 
especially  notions  of  the  2nd  Law  and  entropy,  can  help  clear  a way.  We  ask  if 
such  ideas  offer  practical  means  to  advance  practical  knowledge,  and  where  ma- 
jor impediments  may  lie. 


A discussion 

These  pages  follow  notes  taken  throughout  ‘Aha 
Huliko’a,  and  from  ideas  discussed  in  part  during  the 
‘Aha  discussion  time.  The  issues  were  framed  by 
David  Marshall  and  further  by  Peter  Killworth.  Sur- 
veying approaches  that  have  been  taken  to  providing 
a basis  for  eddy  parameterization,  one  is  daunted  by 
the  length  of  the  list  and  the  tentative  character  of  the 
entries.  What  “should”  eddies  do?  On  the  list  we 
find  (as  examples) 

Eddies  should  flatten  isopycnals. 

Eddies  should  maximally  dissipate  APE  (available 
potential  energy). 

Eddies  should  mix  PV  (potential  vorticity). 

Eddies  should  mix  layers  thickness  along  layers. 

Eddies  should  maximally  dissipate  enstrophy. 

Eddies  should  reduce  MKE  (mean  kinetic  energy). 

Eddies  should  relax  toward  certain  rectified  (“nep- 
tune”)  flows. 

. . . and  so  on. 

“Should”  eddies  do  any  of  this?  While  we  pose 
this  list  in  terms  of  eddy  parameterizations,  another 
‘Aha  topic — stratified  mixing — would  generate  yet 
another  list.  Are  we  dispirited?  To  “help,”  Bill 
Young  inserted  an  estimate  of  the  number  of  active 
degrees  of  freedom  in  the  ocean,  suggesting  1028  “per 
cell”  or  1037  if  one  includes  biology.  (There  was  a 
little  dispute  about  numbers,  but  the  key  message  is 
the  numbers  are  “big” — far,  far  bigger  than  modem 
supercomputers  are  able  to  prognose,  which  are  more 
like  107  tolO9.)  Given  such  circumstances,  Walter 
Munk  asked  if  we  deem  the  situation  hopeless. 


Answering  Walter’s  question  depends  upon  what 
ocean  modelling  seeks  to  do.  If  our  modelling  pro- 
ject amounts  to  trying  to  invent  a steam  engine  from 
molecular  dynamics  simulation  of  water  vapour,  it 
may  well  be  hopeless  in  our  lifetimes.  If  we  would 
invent  a steam  engine  based  on  thermodynamic  func- 
tions, in  part  from  empiricism  and  in  part  from  statis- 
tical physics,  it  may  not  be  so  hopeless.  Importantly 
we  need  to  turn  the  huge  number  of  degrees  of  free- 
dom from  threat  to  opportunity. 

When  we  recognize  that  we’ve  no  ability  nor  prac- 
tical interest  to  know  the  ocean  in  all  its  1037  (or 
whatever)  details,  we  naturally  turn  to  probabilities  of 
oceans.  Mel  Briscoe  asked  if  we  would  predict  evo- 
lution of  probabilities  distributions  or  if  we  limit  at- 
tention to  moments  (expectations)  from  those  prob- 
abilities. Framing  issues  in  terms  of  moments  might 
render  the  task  manageable? 

Lessons  in  our  coffee  cups? 

A difficulty  may  be  in  part  “cultural”  insofar  as 
we,  as  a community,  have  little  orientation  toward 
statistical  physics,  basing  ocean  dynamics  instead  on 
the  classical  mechanics  of  GFD  amended  with  sundry 
by-guess-and-by-golly  mixing  coefficients.  During 
‘Aha  two  other  themes  recurred.  We  were  reminded 
of  the  oft-cited  stirring  cream  into  coffee.  And  we 
were  reminded  of  the  influence  Carl  Eckart  brought, 
seeking  to  base  physical  oceanography  upon  underly- 
ing physics.  Although  we  cannot  invite  Carl’s  direct 
input,  we  might  seek  in  a spirit  after  Carl  to  ask  why 
does  cream  in  coffee  turn  brown.  Here  I only  substi- 
tute my  own  comment.  I should  hope  the  answer  is 
not  because  stirring  causes  enhanced  mixing  (diffu- 
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sion).  I should  hope  the  answer  is  that  internal  inter- 
actions within  the  coffee  cup  transition  the  probabili- 
ties of  cream  and  of  coffee  to  a distribution  with 
higher  entropy.  Practically,  the  useful  representation 
of  this  idea  may  well  be  that  stirring  leads  to  en- 
hanced mixing.  We  see  this  as  a result  following 
from  the  underlying  basis,  after  which  we  might 
quasi-empirically  parameterize  cups,  teaspoons, 
manners  of  agitation,  etc. 

In  the  case  of  stirring  cream  into  coffee,  one  is 
quite  inclined  merely  to  nod  to  the  entropy  discussion 
before  proceeding  directly  to  parameterizing  the  stir- 
ring-mixing. Were  the  topic  of  ocean  eddies  this 
simple,  we  would  hardly  speak  of  so  simple  a matter 
at  ‘Aha  Huliko’a.  But  eddies  are  not  simple.  Then 
stirring-mixing  intuitions,  post-hoc  modified  by  crite- 
ria such  as  listed  at  the  outset  of  this  note,  are  notori- 
ously unreliable.  Was  there  something  more  to  learn 
in  our  coffee  cups? 

Dynamics  of  moments  of  probable  seas 

The  following  story  is  not  yet  clear,  in  part  because 
the  methods  are  so  little  explored.  Here  I make  a 
sketch,  indicating  some  research  directions,  practical 
results,  relations  to  issues  David  and  Peter  framed, 
and  current  challenges. 

First  we  embrace  the  idea  that  oceans  are  known 
only  in  probability.  The  detailed  state  of  the  ocean 
might  be  expressed  in  a state  vector  y,  whose  dimen- 
sion could  well  be  1037,  or  whatever.  We  don’t  know 
y.  We  only  speak  of  elemental  probability  dP  = 
P(y)y  that  the  actual  value  of  y falls  within  a phase 
volume  dy  about  any  given  y.  So  far  the  discussion  is 
aerie-faerie.  What  we  really  would  like  to  access  are 
moments  of  P,  such  as  Y = y dP  with  the  integral 
over  all  y.  Importantly,  the  dimension  of  Y need  no 
longer  be  I037.  We  can  “project  out”  as  much  of  y as 
we  care  not  to  consider,  perhaps  taking  as  only 
“lumped”  (space-time  averages)  over  y.  (My  nota- 
tion could  be  embellished!)  Dimensions  of  Y might 
be  only  107  or  103  or  maybe  only  10.  The  key  con- 
sideration is  that  these  Y are  moments  of  probabili- 
ties. That  distinction  can  get  lost  when,  for  example, 
we  look  at  output  of  GCMs  and  see  maps  of  veloci- 
ties, temperatures,  elevations  or  whatever.  Even 
when  we  admit  that  the  Y are  grid-cell-averaged 
variables,  we  tend  still  to  speak  of  “velocity”  rather 
than  “velocity  moment  of  probability  of,”  for  exam- 
ple. So  what? 


The  “so  what”  gets  us  when  we  write  dynamic 
equations.  Too  easily  we  look  up  equations  for  ve- 
locity or  temperature  or  such  from  textbooks,  the 
only  ambiguity  arising  from  nonlinearities  which, 
averaged  over  space-time  volumes  need  some  closure 
“approximations”.  Were  we  to  ask  instead  for  the 
equation  of  motion  of  the  temperature  moment  of  the 
probable  state,  say,  we  might  (1)  grow  tired  and  (2) 
pause  on  our  way  to  the  textbook.  Until  we  are  clear 
what  are  the  dependent  variables  in  the  problem,  as- 
suming equations  of  motion  is  premature. 

Next  steps  are,  in  part,  familiar.  Linear  terms  in 
equations  of  dyldt  commute  with  expectation  and 
averaging  operators,  so  linear  terms  in  dY/dt  are  “as 
usual”.  Nonlinearities  in  dyldt  can  be  expressed  in 
parts  as  corresponding  nonlinearies  among  compo- 
nents of  Y,  which  again  may  look  familiar  following 
“usual”  Reynolds  averaging.  And  there  is  “more,” 
the  “stuff’  that  connects  the  Y to  all  the  P(y)  which 
we  do  not  know.  dY/dt  = f(Y)+X,  where  “f’  express 
the  “familiar”  equations  from  textbooks  and  “X”  are 
the  unknowns  (of  course). 

Two  routes  to  “X” 

The  question  of  “X”  should  be  seen  in  context  of 
nonequilibrium  statistical  mechanic,  a gloriously  un- 
solved problem.  There  are  two  avenues.  I have 
tended  to  follow  Lars  Onsager,  seeing  in  “X”  the 
generalized  thermodynamic  forcing  X=K*VyS  where 
S = -j  In  (P)  dP  is  entropy,  denotes  the  gradient  of 
entropy  with  respect  to  the  Y,  and  K supplies  the  cou- 
pling with  which  VyS  force  dY/dt.  As  we  don’t  know 
P,  hence  we  don’t  know  S,  or  VyS  and  we  don’t 
know  k,  all  this  looks  like  useless  window  dressing. 
Maybe  not.  If  we  can  determine  some  Y = Y*  for 
which  VvS  is  small  (in  the  sense  much  smaller  than 
VyS  at  the  actual  Y)  we  could  try  to  expand  K*VyS  = 
k»V2yyS*(Y=Y*).  Call  k*V2yyS  = C so  it  doesn’t 
look  so  scary  and  we  have  only  two  problems:  what 
is  Y*  and  what  is  C? 

Y*  is  usually  obtained  by  thinking  about  dY/dt  = 
f(Y)  under  idealized  circumstances,  where  we  sup- 
pose many  excited  degrees  of  freedom  while  omitting 
all  external  forcing  and  internal  dissipation  (here  re- 
garded as  “external”  to  dynamics  of  Y).  Dynamics 
sometimes  are  further  simplified,  e.g.  to  quasi- 
geostrophy,  to  make  calculation  of  Y*  tractable. 


PRINCIPLES  OF  EDDY  PARAMETERIZATION 

Subject  to  integrals  of  the  motion  of  idealized  dY/dt= 
f(Y),  Y*  is  the  Y that  maximizes  S,  i.e.,  V\S  = 0. 

Aside:  This  point  has  confused  onlookers  more 
than  any  other.  The  “theory”  appears  to  be  to  maxi- 
mize entropy.  But  such  a result  would  only  apply  to 
a mathematical  idealization  (an  isolated,  unforced, 
nondissipative  system)  arguably  far  from  Earth’s 
oceans.  What  needs  be  emphasized  is  that  the  whole 
idea  is  to  use  Y*  as  a means  to  access  non-zero  W\S 
in  order  to  complete  the  actual  equations  of  motion  of 
actual  Y.  This  is  not  a “maximum  entropy”  theory  of 
anything. 

Couplings  C remain  to  be  estimated  and,  in  my 
work  to  date,  are  largely  fudged.  (Y*  isn’t  so  great 
either.)  At  this  time  the  point  is  not  to  find  “the  an- 
swer” (don’t  I only  wish!),  but  rather  to  identify  the 
parts  of  the  answer  which  may  yield  to  successive 
efforts.  At  C I encounter  the  same  kinds  of  semi- 
empirical,  largely  fudged,  by-guess-and-by-golly  es- 
timations which  are  characteristic  of  our  ability  to 
represent  oceanic  turbulence. 

Here  let  me  mention  a second  approach  to  “X,”  re- 
cently advanced  by  Joel  Sommeria  and  colleagues. 
The  idea  is  to  find  an  expression  for  overall  produc- 
tion of  entropy,  dS/dt,  which  can  be  maximized  with 
respect  to  Y.  X is  then  the  force  on  dY/dt  which 
maximizes  dS/dt.  Although  both  the  derivation  of 
dS/dt  and  the  assignment  of  constraints  for  maximiz- 
ing dS/dt  have  raised  new  issues  and  new  uncertain- 
ties, the  maximum  entropy  production  approach  of- 
fers an  important  complement  to  the  entropy  gradient 
forcing  which  I have  pursued.  Happily,  Igor  Polya- 
kov has  compared  the  two  approaches  in  a case  of 
Arctic  ocean  modeling  and  finds  pleasingly  similar 
results. 

So  what? 

A reader  can  well  ask:  if  we  are  only  stirring  cream 
into  coffee,  isn’t  this  entropy  talk  a lot  of  bother? 
Indeed  if  only  ocean  eddies  were  as  simple  as  the 
coffee,  we  should  hardly  bother.  For  much  of  what 
we  do  about  ocean  mixing  parameterizations,  effort 
to  recast  the  discussion  in  terms  of  entropy  calculus 
would  (most  likely)  only  append  a superstructure 
over  what — practically — we  do  anyway.  On  the 
other  hand,  during  this  ‘Aha  there  were  two  very  dif- 
ferent areas  of  research  where  results  were  not  “sim- 
ple” in  the  sense  of  coffee  turning  brown. 
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First  recall  David’s  results  in  two-layer  flow  over 
topography.  Numerical  experiments  do  not  lead  to 
flattening  isopycnals  (reducing  APE)  and  do  not  lead 
to  uniform  PV.  What  are  eddies  doing?  The  sugges- 
tion, which  would  need  to  be  quantified  suing  actual 
code  for  actual  geometry  of  David’s  experiments,  is 
that  eddies  move  the  two-layer  flow  to  nearly  the 
highest  entropy  it  can  attain.  One  approach  may  be, 
if  the  code  David  used  can  be  run  in  dissipationless, 
conservative  mode,  then  the  model  itself  can  be  let 
run  to  reveal  Y*.  There  is  no  reason  Y*  should  re- 
flect either  minimum  APE  or  uniform  PV.  When 
actual  dissipation  and  forcing  (if  present)  cause  actual 
Y to  depart  from  Y*,  eddy  fluxes  should  arise  in  the 
model  (testably)  approximately  proportionally  to  Y*- 
Y. 

Peter  reminded  us  of  an  older  illustration  from  sta- 
tistical mechanics,  recalling  a hypothetical  Arctic 
circulation  (from  myself  from  ‘Aha  Huliko’a,  1993!) 
in  which  rectified  (“neptune”)  flows  were  induced  by 
eddies.  While  those  early  results  were  barotropic, 
extensions  to  baroclinic  flow  apply  to  David’s  case. 

A different  result  that  stirred  controversy  during 
‘Aha  was  George  Camevale’s  simulation  of  internal 
wave  breaking.  When  George  evaluated  vertical 
buoyancy  flux,  w'b',  where  buoyancy  b = (p0-p)/p0  is 
the  fractional  deficit  of  density  about  reference  p0  and 
w is  vertical  velocity,  spectral  contributions  were 
positive  (upwards)  over  nearly  all  k.  In  particular 
w'b'(k)  > 0 over  all  k that  were  “turbulent”  by  any 
measure  of  “turbulence”.  Because  the  experiments 
were  stably  stratified,  mean  db/dz  > 0 and  the  turbu- 
lence from  internal  wave  breaking  forced  b up  the  b- 
gradient  (on  average),  “anti”-mixing.  This  is  not  stir- 
ring the  cream  into  the  coffee!  What  was  wrong? 
Sentiments  at  ‘Aha  ranged  from  (1)  the  experiments 
were  performed  improperly  (wrong  large  scale  forc- 
ing) to  (2)  analyzing  outcome  in  “z”  is  wrong,  and 
density  coordinates  should  be  used. 

Or  maybe  George  had  things  right,  as  indeed  (I 
think)  explains  the  differential  diffusion  which  I re- 
port elsewhere  in  these  proceedings.  If  George  was 
right,  why  w'b'(k)  > 0?  Again  I’ll  only  speculate 
without  direct  access  to  George’s  output,  but  I be- 
lieve that  wave  “breaking”  efficiently  scatters  poten- 
tial energy  PE  (as  b'2)  to  higher  wavenumbers.  Over 
most  k,  the  result  was  PE  > KE  at  each  k.  In  these 
experiments  without  Coriolis,  Y*  for  internal  gravity 
modes  equipartitions  PE  and  KE.  Thus  w'b'(k)  > 0, 
converting  PE  =>  KE,  is  driven  by  Y-Y*. 
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Does  this  help?  A shortened  list  of  guiding  princi- 
ples can  read  only  “dS/dt  > 0”  but  the  ongoing  practi- 
cal challenge  is  to  put  this  idea  to  work.  Progress  is 
slow  because  the  methods  are  unfamiliar.  But  I think 
tangible  practical  progress  is  being  made. 


